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ABSTRACT 

We study the atomic neutral hydrogen (HI) content of ~I600 galaxies up to z ~ 0.1 using stacking techniques. The observations were 
carried out with the Westerbork Synthesis Radio Telescope (WSRT) in the area of the SDSS South Galactic Cap (SSGC), where we 
selected a galaxy sample from the SDSS spectroscopic catalog. Multi-wavelength information is provided by SDSS, NVSS, GALEX, 
and WISE. We use the collected information to study HI trends with color, star-forming, and active galactic nuclei (AGN) properties. 
Using NUV - r colors, galaxies are divided into blue cloud, green valley and red sequence galaxies. As expected based on previous 
observations, we detect HI in green valley objects with lower amounts of HI than blue galaxies, while stacking only produces a 3-cr 
upper limit for red galaxies with Mhi < (5 ± 1.5) x 10* Mq and Mhi/Lj < 0.02 ± 0.006 Mq/Lq (averaged over four redshift bins up 
to z ~ 0.1). We find that the HI content is more dependent on NUV - r color, and less on ionization properties, in the sense that 
regardless of the presence of an optical AGN (based on optical ionization line diagnostics), green-valley galaxies always show H I, 
whereas red galaxies only produce an upper limit. This suggests that feedback from optical AGN is not the (main) reason for depleting 
large-scale gas reservoirs. 

Low-level radio continuum emission in our galaxies can stem either from star formation, or from AGN. We use the WISE color-color 
plot to separate these phenomena by dividing the sample into IR late-type and IR early-type galaxies. We find that the radio emission 
in IR late-type galaxies stems from enhanced star formation, and this group is detected in H I. However, IR early-type galaxies lack any 
sign of HI gas and star formation activity, suggesting that radio AGN are likely to be the source of radio emission in this population. 
Future HI surveys will allow for extending our studies to higher redshift, and for testing any possible evolution of the HI content in 
relation to star-forming and AGN properties up to cosmologically significant distances. Such surveys will provide enough data to test 
the effect of radio/optical AGN feedback on the HI content at lower, currently rather unexplored HI detection limit (Mhi < lO’ Mq). 

Key words. Keywords should be given 


1. Introduction 


The amount and conditions of cold HI gas in galaxies are, 
in a direct or indirect way, related to star formation (SF) 
processes and to black hole fuelling. Therefore, our knowl¬ 
edge of the HI properties is crucial for our understanding 
of the intricate processes of galaxy formation and evolu¬ 
tion. Our knowledge of the gas content in various types of 
galaxies in the nearby Universe has increased substantially 
thanks to large single-dish HI surveys such as t he HI Parkes 
All S ky Survey (HIPASS) (iMever et al .L 12004 |z waan et all 


dGiovanelli et al.Ll2005UHavnes et al.. 201 11) and detailed imag- 

ing surveys like WHISP dvan c 

er Hulst. van Albada. & Sancisi 

2001 

), THINGS dWalter et al., 

2008), SAURON dMorganti et al. 

2006 

;IOosterloo et al.ll20I0l) and ATLAS^^ dSerra et al.. 2012). 


Scaling relations linking the H I-to-stellar mass ratio to 
global physical properties of galaxies, e.g. stellar mass, stel¬ 
lar mass surface density and colour, have been derived for 
massive galaxies in t he Galex Arecibo S DSS Survey (GASS, 
ICatinella et al.|[2M^ ICortese et alJl201 ih . These studies show 
that the HI gas fraction most strongly correlates with NUV - 
r color and stellar surface mass density, which relation is 


most naturally explained as being t he Kennicutt-Schmidt law 
(lKennicuttLll998t|Zhang et al.Ll2009l) . The correlation is not sur¬ 
prising, since the NUV - r color and ju* correlate with spe¬ 
cific star formation rate and morphological type respectively, 
which properties are expected to be closely tied to gas frac¬ 
tion. Blue galaxies, in general, are HI rich and actively form¬ 
ing stars, whereas red galaxies tend to be more H I-poor or de¬ 
pleted of gas. The ‘green valley’ is considered to be a transition 
pop ulation between star-forming blue and quiescent red galax¬ 

ies (iMartin et al.L[20071: IWvder et al.L 120071: ISchiminoviche^al . 
2007), displaying residual star-formation signatures (lYietak . 
20051) . As a results, green, intermediate colors have been inter¬ 
preted as a sign for the recent q uenching of star formation in this 
population (ISalim et al.Ll2007h . 


The mechanisms that drive quenching processes in galaxies 
are yet to be understood. It is thought that feedback processes are 
important in affecting the gas reservoirs and consequently the 
star formation processes, particularly in massive, bulgy galax- 
ies ( Crotonl [2006t IPe Lucia & Blaizoll 120071 : l^merville et al.L 
I 2 OO 8 ). Theoretical studies emphasize the role of AGN in 


quenchi n g star formation processes by heatin g the gas rese 
dCrotonl T 2 OO 4 iDe Lucia & Blaizol l2007t ISornerville 


heatin g the gas reservoirs 
' et all 
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l2008t iBooth & Schavei [200^ iDebuhr et all l2Q12h. or by driv¬ 


ing the gas outside of the host gakxy j[pi Matteo et all l2005t 
iHopkins et al.L 120061 : IWagner et all I2OI2I) . Therefore, observa- 
tional evidence is of key importance for our understanding of 
the effects that AGN can have on the (neutral) gas and star for¬ 
mation processes in galaxies. 

Active galactic nuclei (AGN) can reveal their presence in 
many ways, e.g., optical emission lines in the nuclear region, 
and radio jet activity. AGN are often selected optically, based 
on Baldwin, Phillips & Terlevich (BPT) [OIII1/H^ vs [N IIl /Ha 
line ratio d i agnos tics dBaldwin, Phillips. & TerlevichLll98ll) . As 
iBest et ai] (|2005|) pointed out, optical line selection naturally 
leads to a bias in favour of AGN with a rich inter-stellar medium. 
This happens for the reason that to optically identify an AGN, 
the presence of black holes and sufficient amounts of gas that 
produce the bright emission are both necessary. However, AGN 
which are selected in the radio continuum often show only weak 
or no emission lines, and this population would be missed by 
optical selection. 

Radio emission in galaxies can either stem from AGN ac¬ 
tivity, or it can also originate from star formation processes. 
Therefore, care has to be taken to separate these phenom¬ 
ena when looking for AGN in the radio continuum. From the 
1.4 GHz radio luminosity function it is expected that radio 
AGN become dominant over star formation at radio power 
log(P / W Hz“') > 23 (Mauch and Sadler 2007). As a result, the 
highest radio power sources are most likely associated with ac¬ 
tive galactic nuclei (AGN). However, the nature of the faint radio 
population is controversial. Low radio luminosity sources have 
been identihed as being either SF galaxies (starbursts, spirals 
or irregulars), or low power radio-loud and/or radio-quiet AGN 
(e.g. faint FR 1, Seyfert galaxies). Many different approaches 
have been taken to identify SF and AGN processes in low lumi¬ 
nosity radio sources (see Prandoni et al. 2009 for an overview). 
Among these, IR color-color diagrams were found to be efficient 
in disentangling these two phenomena. 

In order to learn about the SF and AGN processes of dif¬ 
ferent galaxy populations in relation to their HI gas content, we 
need large samples to globally look at the HI properties of galax¬ 
ies. Large samples can be achieved either by large sky surveys 
such as ALFALFA, HIPASS, or by looking at smaller portions 
of the sky while extending the observations to higher redshift. 
The limitation of the second approach is that at higher redshift, 
HI studies are limited by sensitivity and bandwith, and such 
observations are quite expensive in terms of observation time 
even with the most sensitive telescopes dCatinella & Cortes3 . 
l2015h . However, alternative techniques, like spectral stacking, 
have been efficiently used to measure the average HI con¬ 
tent of galaxies in the relatively distant Universe. In combina¬ 
tion with multiwavelength data, stacking has also provided a 
powerful tool to study the neutral gas propertie s in different 
galaxy popul a tions o ver the last years dLah et all 1200^ 12009 ; 
Fabello et al.L 201 lallR IVerheiien et aTL l2007t ' Delhaize et al 


2013HGerebet al.Ll2013ll . 


In a previous paper dGereb et alJl2013L i.e. Paper 1.) we car¬ 
ried out spectral stacking on galaxies selected in the Lockman 
Hole (LH) held. We reported the detection of HI gas not only in 
normal SF galaxies, but a lso in low ionization nuclear emission 
regio n (LINER) galaxies (iKauffmann et al.L l2003t iKewlev et al.L 
l2006t [Best & HeckmanL |2012|) . We show in good agreement 
with SAURON and ATLAS^°, that albeit in lower amounts, HI and 
star formation are present not just in typical SF galaxies (gener¬ 
ally blue, late type spirals), but also in galaxies with older stellar 
populations (redder), or in galaxies hosting AGN. In Paper 1. we 


successfully separated red (g - r > 0.7) star-forming galaxies 
from potential low-luminosity radio AGN using IR diagnostics. 
We found that the two groups are well separated in terms of HI 
properties as well. The main limitation of the LH study was that 
our selection of different sub-samples was restricted by the small 
sample size. 

In this paper we present the HI properties of a larger sam- 
ple, ~ 1600 galax ies. We conhrm several trends derived by 
iGereb et aP d2013h and expand on these results. This work 
is made possible by the increase, by more than an order of 
magnitude, in the number of galaxies, allowing us to lower 
the detection limit by a factor of 4 using stacking techniques. 
Furthermore, the relatively high spatial resolution (see Section|2]) 
of the observations used in this study reduces the risk of confu¬ 
sion from companion galaxies, often present in previous, single¬ 
dish stacking experiments. The larger sample makes it possible 
to study in more de tail the HI content o f different galaxy popula¬ 
tions. In addition to lGereb et al.l(l2013l) . we study galaxies where 
quenching and feedback are thought to be affecting the gas reser¬ 
voirs, e.g. green valley objects, AGN (selected based on optical 
emission lines). 

A sub-sample of our galaxies is associated wi th radio emis- 
sion in the NRAO VLA Sky Survey (NVSS, ICondon et al.l 
Il998h . We explore the SF and AGN properties of the radio 
sources in relation to their HI content, and we discuss the possi¬ 
bility of separating these phenomena using IR colors. 

In a second paper, using the same dataset, the Qhi will be in¬ 
vestigated in the same redshift range (Hoppmann et al., in prep.). 

In this paper the standard cosmological model is used, with 
parameters Q,„ = 0.3, A = 0.7 and Hq - 70 km s“' Mpc^'. 

2. Observations and sample selection 

The HI observations were carried out with the Westerbork 
Synthesis Radio Telescope (WSRT) at 1.4 GHz, in the period 
May 2011 - October 2012, in the area of the SDSS South 
Galactic Cap (SSGC). Between the coordinates 21^‘ < RA < 2^\ 
10° < DEC < 16° (J2000), 35 WSRT pointings were observed 
and used for HI spectral stacking. The SSGC was chosen be¬ 
cause in this sky area the oversubscription of the WSRT is usu¬ 
ally the lowest. 

The redshift range 0 < z < 0.12 is covered by 8 x 20 MHz 
bands with 128 frequency channels in each band (1280 - 1420 
MHz, the bands overlap 3 MHz). The corresponding velocity 
resolution is 38 km s '. The integration time is 12 hours for most 
of the observations. The synthesized beam is typically 70 x 9 
arcsec. The elongated beam shape is (mainly) the result of the 
low declination observations with the East-West WSRT array. 
However, the wealth of ancillary data available for this region of 
the sky outweighed the disadvantages of having such an elon¬ 
gated beam. 

We use the Sloan Digital Sky Survey (SDSS, York et al. 
2000) to dehne the central coordinates of the 35 WSRT point¬ 
ings. These were selected in the sightline of high density SDSS 
regions with the goal of maximizing the number of detec¬ 
tions. Because of the relatively large held of view and redshift 
range, both low and high density environments are well sampled. 
The observed spectroscopic sample is cross-correlated with the 
Galaxy Evolution Explorer (GALEX, Martin et al. 2005), pro¬ 
viding 1595 galaxies that can be used for stacking . We also use 
the Wide-Eield Infrared Survey Explorer (WISE, IWright et al.l 
I 2 OIOI) to obtain infrared (IR) data for our galaxies. The WISE 
sample is complete to the 99 percent level. In addition, a few 
(50, within a search radius of 15 arcsec) galaxies are identihed 
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Fig.l. 1. Color-magnit ude diagram ( left) 2 . UV-optical color-color plot (right); The galaxies are color coded according to the color- 
magnitude selection bv lWvder et alJ (l2007h . i.e. blue cloud, green valley and the red sequence. NVSS radio sources are marked by 
yellow squares. 




O 

O) 

o 


2.0 

1.5 

1.0 

0.5 

0.0 

-0.5 

- 1.0 


Raidio 

• SF 

• •. Composite 

• • • Optical AGN 


□ ;• • 

» ‘‘to .. H • 

•/ ■ 

tjLsjl*3fi ff ■ ^ 

* * 

145^7^ . • 

a* ^ • 

V 

' • 

1 


-1.0 -0.5 0.0 

log{[NII]/Ha) 


0.5 



0.0 0.5 1.0 1.5 2. 

WISE [4 


0 2.5 3.0 3.5 4.0 4.5 

.6] - [12] (mag) 


Fig. 2. 1. BPT diagram (left). The dashed line (iKauffmann et al.L l2003h is separating SF galaxies (b elow the dashed line) from 
composite galaxies (between dashed and dotted line). Optical AGN are located above the dotted line (iKewlev et al.L 120011 l2006l) . 
Inactive galaxies do not appear in the diagram. 2. WISE IR color-color plot (right). The sources are color-coded according to the 
BPT selection. The vertical solid line is separating IR early-type ([4.6/tm] - [12pm] < 2) and IR late-type galaxies ([4.6pm] - [12pm] 
> 2). Radio sources are marked by yellow squares. 


with radio counterparts in the NVSS survey. The collected multi¬ 
wavelength data allows us to combine various galaxy parameters 
and investigate the HI properties of different galaxy populations. 

3. Data reduction and H i stacking 

The data were reduc e d usi ng the MIRIAD package 
dSault, Teuben, & Wrighll Il995h . Bad data were flagged 
from the datasets, with extra care for the most prominent RFI in 
the lowest-frequency band (1280 - 1300 MHz). 

The standard way to subtract the continuum is by htting 
a low-order polynomial to the line-free channels. Our datasets 
cover a broad, 140 MHz bandwidth composed by 8 spectral win¬ 
dows. We fit each spectral band separately; however, in cases 
where the HI line is at the edge of the band, instead of interpolat¬ 
ing we need to extrapolate due to lack of continuum at the edge. 


This makes the polynomial fitting unsuccessful and the contin¬ 
uum subtraction more uncertain. These HI lines create a dip in 
the stacked spectra because of the bad continuum subtraction. 

To avoid this effect, first we perform a continuum subtrac¬ 
tion in the wv-plane, using the clean components of each field. 
The clean components were created as a result of deconvolution 
of the continuum images with the dirty beam. This subtraction 
step removes most of the continuum from the HI data cubes. 
However, we perform a second continuum subtraction by fitting 
a second-order polynomial to the spectra to subtract low-level 
residual continuum emission coming from very bright sources. 
Following these steps, we eliminate the dip from the stacked 
spectra. 

Sta cking is done similarly as described in iGereb et^ 
(1201 3h . centred on the redshift of the galaxy to be stacked. 
Galaxies are stacked in four redshift ranges, between 0.02 < z 
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< 0.12 each redshift range covers Az — 0.02 (except the highest 
redshift range, 0.08 < z < 0.12). After noise corrections in every 
channel, the weighted sum of the stacked spectra is; 


5(y) = 


y S i(y) 

micrjiv)) 


(1) 


where 5(v) is the stacked spectrum, cr,(v) and 5,(v) are the 
channel noise and the extracted spectrum of source i as a 
function of frequency. 


The HI masses are derived from the formula: 

Mhi _ 2.356 X lOV ^yW Dl W AV \ 
Mq (1 +z) \Jy j\Mpc/ \km s“* j 



where z is the mean redshift of the stacked sample, 5v is the 
mean flux integrated over the AV velocity width of the HI profile 
in the emitter’s frame, and Di is the average luminosity distance. 
The r-band is chosen to calculate luminosities for the galaxies, 
which are then used to evaluate Hl-mass to stellar luminosity 
ratios. 

The advantage of stacking the HI flux in small redshift bins 
is that one can track any possible evolution of the HI content as 
function of redshift. However, when exploring the HI properties 
over the entire redshift range, we stack gas fractions instead of 
HI flux. This way we compensate for the fact that, for the same 
amount of H I, the flux emitted by nearby HI sources is stronger. 
The extracted spectrum of source i is expressed in terms of HI 
mass-luminosity fraction: 


5,(y)[/y] 


5;(y)[JyMpc2L-'] 


2.356 X 10^ 


1 +Zi 




(3) 


where zi is the redshift, and Dm is the luminosity distance, and L, 
is the r-band luminosity of source i. After stacking, the mean HI 
mass-luminosity ratio is integrated over the AV velocity width. 
Before stacking, the average nns noise of the cubes is ~ 0.2 
mJy/beam, which is expected to decrease with the square root 
of the number of co-added sources. 


4. Characteristics of the gaiaxies in the seiected 
sampie 

Our main goal is to study/compare the HI properties of various 
populations of galaxies using large samples with available mul¬ 
tiwavelength information. In order to do this, we define several 
sub-samples, using the collected multiwavelength data. 

The color distribution of galaxies is known to consist of 
two main peaks, i.e. the blue cloud and the red sequence 
(IStrateva et alll20()TtlBaldrv et al.Ll2004l) . However, at fixed ab¬ 
solute magnitude the two peaks are not well fit by a double 
Gaussian due to an excess of gala xies at intermediate colors, 
also known as the ‘green valley’ dWvder et al.L l2007h . Such 
intermediate, green colors can be due to a number of differ- 


ent 
ity ( 

phenomena, e. 

g. low level (residual) star formation activ- 

Yi et a 

. 2005 

), dust extinction, older stellar populations 

(ISarzi et al. 

2010h. 


Previous stacking studies, including our LH analysis, were 
carried out on galaxies separated into blue/red samples, which 
colors were defined based on ultraviolet-optica l NUV - r, 
or optical g - r selections (iFabello et al.L 1201 lat iGereb et"^ 


Fig. 3. Radio power histogram of NVSS sources in our sample. 


l2013h . To expand on previous studies, in this paper we con- 
sider the green v alley as a separate population. Following 
IWvder et alJ (I2nn7h . in Fig. [T] our objects are divided into blue 
cloud, green valley and red sequence objects based on NUV 
- r colors. To der ive the color distri bution o f the g alaxies, 
IWvder et al.l (l2007h utilizes the fit by lYi et al.l (l2005l) to the 
NUV - r colors in function of the Mr absolute magnitude: 
NUV - r = f(Mr) = 1.73 - 0.17M,. The red sequence is de¬ 
fined as the galaxies with NUV - r > f{Mr) - 0.5, blue galax¬ 
ies have NUV - r < f{Mr) - 2, whereas green valley galax¬ 
ies are the excess population between blue and red galax¬ 
ies with f{Mr) - 2 < NUV - r < f{Mr) - 0.5 colors. The opti¬ 
cal and UV apparent magnitudes ar e extracted from GALEX 
and SPSS, and K-corrected follow ing IChilingarian et all (l2010l) : 
IChilingarian & ZolotukhinI d2012h . The NUV - r colors are cor¬ 
rected for Galactic extinction following lWvder et al.l (l20 07h. 

Fr om Fig. [T] (right) it is clear that in case of the iGereb et al.l 
(12013h selection, our current green valley objects are part of the 
g - r > 0.7 (red) sample, and in the selection of IFabello et al.l 
(1201 lah . green valley objects belong to the NUV - r < 4.5 (blue) 
sample. 

We also use the Bald win, Phillips & Terlevich (BPT) line 
ratio diagnostic diagram dBaldwin, Phillips, & Terlevichl Il981l) 
to separate galaxies with different ionization properties (Fig. |2] 
left). The selection is done b y using line fluxes from SDSS, 
similar to IGereb et'aP (12013h . but with o ne major difference . 
Galaxies, which were defined as LINERs in IGereb et alJ(l2013ll . 
are now separated into composite galaxies and optical AGN 
sampl es, using the more stringent demarcation bv iKewlev et al.l 
(l2006h for selecting AGN. Following the BPT classification, our 
sample includes star-forming (SF) galaxies, composite galax¬ 
ies where the optical line excitation can come from both star 
formation and AGN activity, and optical AGN consisting of a 
mix of Seyfert and Low Ionization Nuc l ear Emission Region 
(LINE R) galaxies (iKewlev et al.L l2006t ICid Eernandes et ^ 
120101) . Eurthermore, a part (280 sources) of our sample is de¬ 
fined as optically inactive, with non-detected, or with maximum 
two detected lines (among [Nil], Her, [OIII], H/?). Star-forming 
galaxies in the BPT diagram in Pig. |2] (left) are mainly blue and 
green, while composite galaxies and optical AGN are green and 
red. Inactive galaxies are typically red, as in Paper 1. 

The NVSS radio population is marked by yellow squares in 
all panels of Pig. [1] Radio sources are among the optically most 
luminous objects in each population of the color-magnitude di- 
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Color 

Redshift bin 

Mhi (10^ Mq) 

Mni/Lr (Mq/Lq) 


0.02 < z < 0.04 

2.22 ± 0.05 

0.45 ± 0.009 

g - r <0.7 

0.04 < z < 0.06 

4.04 ± 0.07 

0.48 ± 0.008 

0.06 < z < 0.08 

4.71 ±0.12 

0.38 ± 0.009 


0.08 <z < 0.12 

7.49 ± 0.24 

0.39 ±0.012 


0.02 < z < 0.04 

0.76 ± 0.04 

0.06 ± 0.003 

g-r >0.7 

0.04 < z < 0.06 

1.42 ±0.07 

0.08 ± 0.004 

0.06 < z < 0.08 

1.81 ±0.09 

0.07 ± 0.004 


0.08 <z < 0.12 

3.01 ±0.21 

0.10 ±0.007 


Table l.Hl mass and mass-luminosity ratio derived for blue (g - r < 0.7) and red (g - r < 0.7) galaxies in the four redshift bins. 




z 


Fig.4. (Left): HI mass-luminosity ratio of galaxies in our current color selection of blue cloud, green valley and red sequence 
objects. Red arrows and dashed lines indicate the 3-cr upper limits. (Right): The SFE of blue cloud/green valley objects. The blue 
dashed line indicates the region of high efficiency star formation, whereas the low efficiency region is marked by a red dashed line. 
The error on the SFR is estimated from the 3-(T NUV flux errors. 


agram in Fig. [T] (left). From the radio power distribution in Fig. 
[3] we expect to have powerful AGN at log(P / W Hz *) > 23, 
however the distribution shows that low-power radio sources are 
also present in the sample. In the low radio power regime it be¬ 
comes more complicated to disentangle the contribution of SF 
and AGN to the radio continuum emission. 

As previous multiwavelength analysis has shown, IR colors 
are efficient in sep arating SF and AGN processes in galaxies 
(iGereb et al.ll2013l and references therein). With the goal of dis¬ 
entangling SF and AGN activity in our radio sample, we extract 
3.4 fim, 4.6 jum and 12 /urn magnitudes from WISF to constrain 
the IR color-color plot, presented in Fig. |2] (right). The sepa¬ 
ration at the vertical line (at [4.6/tm] - [llfim] = 2) in the IR 
color-color plot (Fig. |2] right) is often use d in the literature to 
disentangle IR earl y- and late-type galaxies (I Wright et al.Ll2010l: 
ISadler et al.L 1201 3h . or in other words non-star-forming and SF 
galaxies. Correspondingly, the radio sources are separated into 
two distinct populations in Fig.|2](right). A more detailed analy¬ 
sis of this separation is presented in Sec. 15.31 

5. Results 

5.1. Stacking in coior 

Before we look at the HI content of blue cloud/green val¬ 
ley/red sequence objects, first we compare our current mea¬ 
surements with the results from the LH study. In the redshift 
range 0.06 < z < 0.09, in the previous LH held we mea¬ 
sured Mhi = (6.1 + 0.4) X 10‘^ Mo and Mhi/L, = 0.38 + 0.02 


Mq/Lq in the blue (g - r < 0.7) population, whereas red 
{g - r > 0.7) galaxies contain lower amounts of gas, with 
Mhi = (1-8 + 0.2) x 10® Mq andMni/L.. = 0.08 + 0.01 Mq/Lq. 

In Table [T] we use g — r optical colors to evaluate the HI 
mass-luminosity ratio of blue/red galaxies in our current, larger 
sample. We can average the HI results over the last two red¬ 
shift bins (0.06 < z < 0.12) to achieve a similar redshift selec¬ 
tion as for the LH. In this redshift range, blue galaxies show 
Mhi = (6.1 + 0.2) x 10® Mq and Mhi/L^ = 0.38 + 0.01 Mq/Lq, 
whereas red galaxies have Mhi = (2.4 + 0.15) x 10® Mq and 
Mhi/Li = 0.08 + 0.01 Mq/Lq. These new measurements con¬ 
firm the results that we obtained in the LH study. 

We now move to t he more detailed grou ping we can do with 
the expanded sample. iFabello et al.l (1201 lal) reported HI detec¬ 
tion in red galaxies with NUV - r > 4.5 (similar to our selec¬ 
tion of red galaxies), albeit in lower amounts compared to blue 
galaxies. In Fig. |4] (left), we present the HI mass-luminosity ra¬ 
tio for blue cloud, green valley and red sequence objects sep¬ 
arated based on NUV - r color in our sample. We detect HI 
in blue and green objects, however unlike in the IFabello et al.l 
(1201 lah study where more galaxies are available for stacking, 
red galaxies do not show an HI detectiorQ at low 3-cr limit of 
Mhi < (5 + 1.5) X 10*^ Mq and Mm/Lr < 0.02 + 0.006 Mq/Lq 
( with values averaged over the four redshift bins). As we ex¬ 
pected based on previous studies, we And that green valley ob- 

* We note that in the first redshift hin we find a tentative HI detection 
in red galaxies at the 3-cr level, with Mhi < (3.8 ±0.6) x 10® Mq and 
Mhi/Li < 0.03 ± 0.005 Mq/Lq. Higher redshift bins are not detected. 


5 











































Gereb et al.: The global HI content of SF galaxies and AGN 
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Fig. 5. (Left): The [OIII] luminosity distribution of composite galaxies and optical AGN with green and red colors. (Right): HI 
mass-luminosity ratio of green/red composite galaxies and optical AGN. The 3-cr upper limits are marked by a dashed line and 
arrows. 


jects are a transition population from HI point of view, showing 
lower amounts of HI than the blue population, however they are 
more H I-rich than the red sample. 


Considering that green valley galaxies are expected to have 
lower amounts of HI than the blue population, it is interesting to 
explore the timescale over which green valley objects would de¬ 
plete their gas reservoirs due to star-forming processes. Naively, 
one could expect that the short transition phase of green val¬ 
ley objects is due to quick consumption of their small gas reser¬ 
voirs. However, recent results show that in the nearby Universe, 
the H I-based star formation efficiency (SFE = SFR/Mhi) - or 
the equivalent inverse, the time scale of neutral gas consump¬ 
tion (t = Mhi/SFR) - is independent of other galaxy properties, 
such as stellar mass, stellar surface density, co lor, concentration 
dSchiminovich et ^bOlOUBigiel et~^l2()l ih . 


Here we test the efficiency of star formation in blue/green 
galaxies in our sample, baring in mind their HI content. The star 
formation rate (SFR) is de rived from the NUV flux, following 
ISchiminovich et aP ( 2010t) . The SFR formula accounts for dust 
attenuation by combining UV-optical colors (NUV - r) and the 
D„(4000) index of galaxies. The latter index is an indicator of 
the presence of young stellar populations. In each redshift bin, 
the H I-based star formation efficiency is defined as the average 
SFR over stacked HI mass, i.e. <SFR>/<Mhi>. 


In Fig. |4] (right), green valley objects and blue galaxies re¬ 
veal similar, efficient star formation, with little variations around 
SFE = 10'®'^ yr *, corresponding to a gas consumption time 
scale of t ~ 3x10^ yr. Red galaxies (not shown) lack any sign 
of HI gas and star formation activity. Our estimates for the 
SFR and SFE agree w e ll wit h the results of previous stud¬ 
ies (ISchiminovich et ^ l2010l) . The similar SFE detected in 
blue/green galaxies suggests that the H I-based efficiency of star 
formation is independent of color. 


Our current study spans quite a large redshift range with re¬ 
spect to the majority of previous HI studies. Stacking our galax¬ 
ies in redshift bins allows us to trace the global HI properties as 
function of redshift. We note that the global HI content (mass- 
luminosity ratio) does not change significantly up to z ~ 0 . 1 . 
Even though the sampled redshift range here is too small to relate 
our results to the cosmological evolution of H I, the constancy of 
the global HI content as function of redshift up toz ~ 0.1 is con- 
sistent with what previous studies have found dFreudling et al.L 


1201 ll: iDelhaize et al.L l2013l) . Similarly to the behaviour of the 
global HI content, the SFE shows only little variations with red¬ 
shift up to z ~ 0.1. Over the probed redshift range, the gas con¬ 
sumption timescale of ~1 - 3 Gyr is comparable to the lookback 
time of ~ 1.5 Gyr. Given that the timescales are very similar, gas 
consumption by star-forming processes should have a detectable 
effect on the global HI content. However, the constancy of the 
HI mass-luminosity ratio and SFE as function of redshift is sug¬ 
gestive that the gas in galaxies is replenished on a timescale of 
a few Gyr-s. Similar conclusions for the need of gas replenish¬ 
ment have been reached fo r young galaxies between z = 2 and 
z = 1 (iTacconi et al.Ll2010l) . 

5.2. The HI properties of composite gaiaxies and opticai 
AGN 

Here we want to test the effects that are responsible for the de¬ 
creased HI content in optically selected AGN. Is the lower HI 
content in composite galaxies/AGN the result of galaxies being 
redder (with older stellar populations, dustier), or it is related to 
ionization/AGN feedback properties? 

In the following analysis we select a combined sample of 
composite galaxies and optical AGN considering all galaxies 
above the dashed line in the BPT diagram in Fig. (left), 
and we use NUV - r colors to test the color dependence 
of HI in the selected populations. The strength of the [OIII] 
emis sion line is a fai r ly rel iable tracer of black hole accre¬ 
tion (iKauffmann et al.L l2003ll . To achieve a uniform emission 
line selection, we consider galaxies with [O III] luminosities be¬ 
tween 5 < log(L[oiii]/L 0 ) < 6.5. Following this selection, less 
than 10% of the AGN are excluded from the sample. We exam¬ 
ine the color of the selected objects, and we find that they are 
mostly green and red. In Fig. |5] (left) we show the [OIII] lumi¬ 
nosity distribution of the selected green/red AGN as function of 
redshift. 

To test the HI properties of our AGN as function of color 
and nuclear activity, we stack the gas fractions of red and green 
AGN separately as function of their [OIII] luminosity. The 
stacking is performed in three luminosity bins with steps of 
Alog(L[oiii]/LG) = 0.5. The results of this stacking experiment 
are shown in Fig.|5](right). The HI is concentrated in green ob¬ 
jects among composite galaxies/AGN, whereas red galaxies are 


6 









































Gereb et al.: The global HI content of SF galaxies and AGN 


not detected at the 3-cr level of Mhi/Lj < 0.03 + 0.01 Mq/Lq. 
This suggests that even those galaxies in which the presence of 
an AGN is expected to be more likely (composite galaxies/AGN) 
do show HI detections, however this depends on their color. 
Furthermore, the HI remains relatively constant as function of 
[O III] luminosity. This is a strong indication that the HI content 
is well correlated with the NUV - r color and the SF history of 
the galaxies, while the effect of optical AGN feedback on the dif¬ 
fuse atomic gas is less signi ficant. Similar conclu sion has been 
reache d in previous works bv iFabello et al.l (l201 Ibh and lHo et al.l 

(I2008h . 

5.3. AGN and SF properties of the radio popuiation 

The two main processes that produce radio emission in galaxies 
are AGN and SF. In the nearby Universe, the largest black holes 
are hosted by massive galaxies, which generally have older stel¬ 
lar populations and are poor in gas. Optically these are not very 
powerful AGN, which is also reflected by the fact that about 40% 
of our radio AGN sample is optically inactive with [O III] emis¬ 
sion lines below the 3-cr significance level. However, radio se¬ 
lection is sensitive to detecting this type of AGN. Because radio- 
loud AGN can not be identified based on emission lines, we need 
to look for other multiwavelength diagnostics. 

As we show in Sec 15.21 radio sources are divided into two 
well separated groups in the IR color-color diagram, e.g. IR late- 
type galaxies (on the right) and IR early-type galaxies (on the 
left). In Fig. |2] (right), the IR late-type region is populated by 
star-forming and composite galaxies, whereas the IR early-type 
region is dominated by optical AGN and non-star-forming (opti¬ 
cally inactive) galaxies. Radio AGN are most likely to be found 
in the latter sample. 

As radio AGN are expected to be more powerful than SF 
galaxies, in Fig. |6] we compare the radio power distribution of 
IR late-type and IR early-type radio sources. The two distri¬ 
butions have different shapes, with a wide, D = 0.504 max¬ 
imum distance between the two cumulative distribution func¬ 
tions. According to the Kolmogorov-Smirnov test, the proba¬ 
bility that the two distributions are different is 99%, implying 
that statistically IR late-type and IR early-type galaxies have 
a different radio power distribution. The mean radio power of 
IR late-type galaxies is log(P / W Hz“') = 22.5, whereas as ex¬ 
pected, IR early-type galaxies are typically more powerful, with 
a mean log(P / W Hz“') = 23. 

We want to test whether the radio AGN candidates and opti¬ 
cally selected AGN population show similar properties in terms 
of [O III] luminosity and black hole accretion, or the two AGN 
phenomena are not related. In Fig.|7]we compare the [GUI] lu¬ 
minosity distribution of the two samples using the Kolmogorov- 
Smirnov test. For the AGN which are optically inactive, the 3-cr 
upper limit is considered when calculating the cumulative dis¬ 
tribution. The probability that the two samples are d rawn from 
the sa me distributions is 0.01%. In agreement with iBest et alJ 
(I2OO5I) . this result implies that radio AGN luminosity and opti¬ 
cal AGN emission lines are independent from each other. 

In Paper 1. we carried out HI stacking on radio sources 
separated into SF and AGN samples in a similar way, by IR 
colors. We found that radio-emitting SF galaxies contain HI, 
whereas we did not detect HI in low-luminosity radio AGN at 
the Mhi < (1.2 + 0.40) x 10^ Mq upper limit. Now we do the 
same analysis with radio-selected NVSS objects. 

Stacking results in an HI detection in IR late-type 
galaxies, showing Mhi = (3.04 + 0.37) x 10® Mq and 
Mhi/L, = 0.11 + 0.01 Mq/Lq (the values are averaged over the 



Fig. 6. Cumulative fraction of the radio power distribution in 
IR late-type and IR early-type galaxies. D is the maximum dis¬ 
tance between the cumulative fraction of the two distributions. 
We measure D = 0.504 for 35 IR late-type, and 15 IR early-type 
galaxies. 



Fig. 7. Cumulative fraction of the [O III] luminosity distribution 
of all optically selected AGN (red and green composite-!-AGN) 
in black, and radio-selected AGN in blue. Red and green optical 
AGN are also marked separately. We measure the maximum dis¬ 
tance between the distribution of 321 optical AGN and 14 radio 
sources as D = 0.564. 

entire redshift range). The HI mass-luminosity ratios of these, 
likely star forming galaxies is lower with respect to the en¬ 
tire blue SF population in Fig. |4] The low mass-luminosity ra¬ 
tio is partly the result of the high optical luminosity of radio- 
detected objects, which are among the brightest sources in the 
color-magnitude diagram in Fig. [T] (left plot). It is interesting to 
note that the IR late-type radio population reveals a high SFE 
of ~10 ® yr'*, corresponding to gas consumption time scales 
of t = 10® yr. This SFE is more enhanced than what we detect 
for the entire SE population. Likely, the radio emission in these 
sources is the result of their enhanced star formation activity. 

In the IR early-type region stacking reveals an HI non¬ 
detection with a mass and mass-luminosity upper limit of 
Mhi < (1-95 + 0.6) x 10® Mq andMHi/U- < 0.03 + 0.01 Mq/Lq. 
Along with the high radio power and lack of star formation 
(based on their red colors and location in the IR color-color dia- 
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gram), the HI non-detection supports that radio emission in the 

IR early-type region is due to radio AGN. 

6. Discussion and summary 

• Our stacking results show that galaxies in the green valley 
are detected with lower amounts of HI than blue galaxies, 
but unlike red galaxies, they are not completely depleted of 
cold (HI) gas. This result can be expected based on earlie r 
studies of HI scaling relations (ICatinella et al.Ll^lOll2012l) . 
This result suggests that SF goes hand in hand with the HI 
properties, and in galaxies where cold HI gas is present, con¬ 
ditions are favourable for (residual) SF to be seen. 

• In Fig. 0] the HI mass-luminosity ratio does not change sig¬ 
nificantly as function of redshift, suggesting that the HI con¬ 
tent remains relatively constant up to z ~ 0.1. Over our rela¬ 
tively small redshift range, this result agrees with other stud¬ 
ies of the global HI content. Furthermore, the SFE displays 
a similar behaviour, remaining relatively constant in the cov¬ 
ered redshift range. It seems that the H I-based SFE does not 
depend on color, meaning that green galaxies consume their 
HI on the same timescale as blue galaxies. The fact that the 
H I-based SFE does not seem to depend on galaxy proper¬ 
ties has been interpreted as an indications that the HI con¬ 
tent and SE are regulated by the same process, e.g. feed back 
effects, galaxy environment (ISchiminovich et akl 201(]ll . On 
the other hand, it has been shown that the H2-based SFE 
does depend on several galaxy properties. Systematic stud¬ 
ies of the atomic and molecular gas content of galaxies show 
that there is a significant scatter between the amount of HI 
and H2, which explains while there are clear trends between 
Mh, and SFR, but far weake r ones between Mhi and SER 
(ISaintonge et al.L 1201 ll l2012h . As such, these studies argue 
that the lack of dependencies of the H I-based SEE may 
rather reflect that there is a wide range in the structural prop¬ 
erties of the H I-rich outer disks of galaxies. Perhaps, the HI 
and SE properties of galaxies depend both on environmen¬ 
tal/feedback effects and structural properties of the gaseous 
disks. Euture, high-resolution HI surveys will make it possi¬ 
ble to study the relation of HI to SF properties as function 
of structural/environmental/feedback effects. Future surveys 
will also allow for tracing the evolution of the global HI con¬ 
tent as function of redshift, by extending HI studies to cos- 
mologically significant distances. 

• We detect HI gas in green optical AGN, suggesting that even 
galaxies with higher ionization properties (composite galax¬ 
ies and optical AGN) do contain neutral gas. However, red 
AGN in our sample are not detected in HI. The HI mass- 
luminosity ratio remains relatively constant as function of 
[GUI] luminosity, suggesting that optical AGN are not the 
(main) reason for depleting gas reservoirs. In agreement with 
previous studies, our results show that the presence of HI is 
better correlated with NUV - r color rather than with ion¬ 
ization properties. 

• In radio sources located in the IR late-type region, small 
amounts of gas are associated with very efficient star for¬ 
mation. The lack of HI and the high level of residual star 
formation suggest that these galaxies recently went through 
an intense star-formation period, and this led to a significant 
depletion of HI in these galaxies. 

• We do not detect any HI gas in radio sources located in 
the IR early-type region ([4.6 fim] - [12 jum] < 2) in the 
WISE color-color plot. The lack of HI gas along with the 
non-star-forming properties and high average radio power of 


log(P / W Hz“*) = 23 suggest that the radio emission in this 
population can not originate from star formation. Therefore, 
radio AGN are likely to be responsible for the radio contin¬ 
uum emission in IR early-type galaxies. 

• Radio-selected AGN have a different [O III] lumi nosity dis¬ 
tributi on than optically defined AGN. Eollowing [Best et aH 
(l2005l) . optical and radio AGN represent two different phe¬ 
nomena, triggered in a different way. Unlike optically se¬ 
lected AGN, the radio AGN population is expected to be 
able to deplete/heat cold gas reservoirs (see Introduction). 
On the red sequence, both the optical and radio AGN pop¬ 
ulations have an HI non-detection, and our current sam¬ 
ple is too small to test any differences in the HI proper¬ 
ties of the two AGN populations at low detection limit. 
However, it will be interesting to explore this at lower 
limit with future surveys. In our current selection, we 
reach a detection limit of Mhi < (5 ± 1.5) x 10* Mq and 
Mni/Lr < 0.02 + 0.006 Mq/Lq (averaged over four redshift 
bins) in the red sample. Even though this is a relatively low 
detection limit, lower HI masses have been detected before 
by direct observations of the SAURON and ATLAS*° sam- 
ples, and by earlier st acking studies of ALEALEA galaxies 
(lEabello et al.L 1201 lallbh . Stacking is a promising technique 
to lower the detection limit and study the relatively unex¬ 
plored <10^ Mq Hi mass regime of galaxies using large 
samples of galaxies. This will be made possible by future 
HI surveys with the next generati on of radio telescopes, 
e.g. Apertif (lOosterloo et akLUOlObh l. the Australian Squar e 
Kilometre Array Path finder (ASKAP, iDeBoer et al.ll200^ . 
MeerKat ([Booth et al.L l2009ll . and the Karl G. Jansky Very 
Large Array (JVLA, e.g. CHILES survey, lEernandez et al.l 
l2015h . 
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